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This paper is devoted to the introduction of a new class of consistent estima- 
tors of the fractal dimension of locally self-similar Gaussian processes. These 
estimators are based on convex combinations of sample quantiles of discrete 
variations of a sample path over a discrete grid of the interval [0, 1]. We derive 
the almost sure convergence and the asymptotic normality for these estima- 
tors. The key-ingredient is a Bahadur representation for sample quantiles of 
non-linear functions of Gaussians sequences with correlation function decreas- 
ing as k~ a L{k) for some a > and some slowly varying function L(-). 



1 Introduction 

Many naturally occuring phenomena can be effectively modelled using self-similar pro- 
cesses. Among the simplest models, one can consider the fractional Brownian motion 
introduced in the statistics community by Mandelbrot and Van Ness (1968). Frac- 
tional Brownian motion can be defined as the only centered Gaussian process, denoted 
by (X(t))t£R, with stationary increments and with variance function v(-), given by 
v(t) = cr 2 \t\ 2H , for all t £ 1 The fractional Brownian motion is an //-self-similar 
process, that is for all c > 0, (X(ci)) tgR = c H (X(t)) teR (where = means equal in 
finite-dimensional distributions) with autocovariance function behaving like 0{\k\ 2H ~ 2 ) 
as |fe| —* +oo. So the discretized increments of the fractional Brownian motion (called the 
fractional Gaussian noise) constitute a short-range dependent process, when H < 1/2, 
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and a long-range dependent process, when H > 1/2. The index H also characterizes the 
path regularity since the fractal dimension of the fractional Brownian motion is equal to 
D = 2 — H. According to the context (long-range dependent processes, self-similar pro- 
cesses,. . . ), a very large variety of estimators of the parameter H has been investigated. 
The reader is referred to Beran (1994), Coeurjolly (2000) or Bardet et al. (2003) for 
an overview of this problem. Among the most often used estimators we have: methods 
based on the variogram, on the log— periodogram e.g. Geweke and Porter-Hudak (1983) 
in the context of long-range dependent processes, maximum likelihood estimator (and 
Whittle estimator) when the model is parametric e.g. fractional Gaussian noise, meth- 
ods based on the wavelet decomposition e.g. Flandrin (1992) or Stoev et al. (2006) and 
the references therein, and on discrete filtering studied by Kent and Wood (1997), Istas 
and Lang (1997) and Coeurjolly (2001). We are mainly interested in the last one, which 
has several similarities with the wavelet decomposition method. Following Constantine 
and Hall (1994), Kent and Wood (1997), Istas and Lang (1997), in the case when the 
process is observed at times i/n for i = 1, . . . , n, this method is adapted to a larger class 
than the fractional Brownian motion, namely the class of centered Gaussian processes 
with stationary increments that are locally self-similar (at zero). A process (X(t))t£R is 
said to be locally self-similar (at zero) if its variance function, denoted by v(-), satisfies 

v(t) = B(X(t) 2 ) = a 2 \t\ 2H (1 + r{t)) , with r(t) = o (1) as \t\ -► 0, (1) 

for some < H < 1. An estimator of H is derived by using the stationarity of the 
increments and the local behavior of the variance function. When observing the process 
at regular subdivisions, the stationarity of the increments is crucial since the method 
based on discrete filtering (and the one based on the wavelet decomposition) essentially 
uses the fact that the variance of the increments can be estimated by the sample moment 
of order 2. We do not believe that this framework could be valid for the estimation of 
the Hurst exponent of Riemann-Liouville's process, e.g. Albs et al. (1999) which is an 
iJ-self-similar centered Gaussian process but with increments satisfying only some kind 
of local stationarity, see Remark 2 for more details. 

Let us be more specific on the construction of the wavelet decomposition method, 
see e.g. Flandrin (1992): the authors noticed that the variance of the wavelet coefficient 
at a scale say j behaves like 2^ 2H ~ 1 \ An estimator of H is then derived by regressing 
the logarithm of sample moment of order 2 at each scale against log(j) for various scales. 
This procedure exhibits good properties since it is also proved that the more vanishing 
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moments the wavelet has the observations are more decorrelated. And so asymptotic 
results are quite easy to obtain. However, Stoev et al. (2006) illustrate the fact that this 
kind of estimator is very sensitive to additive outliers and to non-stationary artefacts. 
Therefore, they mainly propose to replace at each scale, the sample moment of order 2, 
by the sample median of the squared coefficients. This procedure, for which the authors 
assert that no theoretical result is available, is clearly more robust. 

The main objective of this paper is to extend the procedure proposed by Stoev et 
al. (2006) by deriving semi-parametric estimators of the parameter H, using discrete 
filtering methods, for the class of processes defined by (1). The procedure is extended in 
the sense that we consider either convex combinations of sample quantiles or trimmed- 
means. Moreover, we provide convergence results. The key-ingredient is a Bahadur 
representation of sample quantiles obtained in a certain dependence framework. Let Y = 
(Y(l), . . . , Y(n)) be a vector of n i.i.d. random variables with cumulative distribution 
function F, as well denote by £(p) and £ (p) the quantile respectively the sample quantile 
of order p. By assuming that F'(£(p)) > and F"(£(p)) exists, Bahadur proved that as 
n — ► +co, 

with r n = O a , s , (n~ 3 / 4 log(n) 3 / 4 ) . Using a law of iterated logarithm's type result, Kiefer 
obtained the exact rate n -3 / 4 log log(n) 3 / 4 . Extensions of the above results to dependent 
random variables have been pursued in Sen and Ghosh (1972) for <j>— mixing variables, 
in Yoshihara (1995) for strongly mixing variables, and recently in Wu (2005) for short- 
range and long-range dependent linear processes, following works of Hesse (1990) and Ho 
and Hsing (1996). Our contribution is to provide a Bahadur representation for sample 
quantiles in another context that is for non-linear functions of Gaussian processes with 
correlation function decreasing as k~~ a L(k) for some a > and some slowly varying 
function L(-). The bounds for r n are obtained under the same assumption as those used 
by Bahadur (1966). 

The paper is organized as follows. In Section 2, we give some basic notations and 
some background on discrete filtering. In Section 3, we derive semi-parametric estimators 
of the parameter H, when a single sample path of a process defined by (1) is observed 
over a discrete grid of the interval [0, 1]. Section 4 presents the main results: Bahadur 
representations and asymptotic results for our estimators. In Section 5 are presented 
some numerical computations to compare the theoretical asymptotic variance of our 



Hurst exponent estimation using sample quantiles 



4 



estimators and a simulation study is also given. In particular, we illustrate the relative 
efficiency with respect to Whittle estimator and the fact that such estimators are more 
robust than classical ones. Finally, proofs of differents results are presented in Section 6. 

2 Some notations and some background on discrete filter- 
ing 

Given some random variable Y, Fy(-) denotes the cumulative distribution function of 

Y and £y(p) the quantile of order p, < p < 1. If Fy(-) is absolutely continuous 
with respect to Lebesgue measure, the probability density function is denoted by /y(-). 
The cumulative distribution (resp. probability density) function of a standard Gaus- 
sian variable is denoted by <&(•) (resp. (/>(■))■ Based on the observation of a vector 

Y = (Y( 1 ),..., Y(n)) of n random variables distributed as Y, the sample cumulative 
distribution function and the sample quantile of order p are respectively denoted by 
Fy (•; Y) and £y (p; V) or simply by F(-; Y) and £ (p; Y). Finally, for some measurable 
function g(-), we denote by g(Y) the vector of length n with real components g(Y(i)), 
for i = 1, . . . , n. 

A sequence of real numbers u n is said to be O (v n ) (resp. o (v n )) for an other sequence 
of real numbers v n , if u n /v n is bounded (resp. converges to as n — > +oo). A sequence 
of random variables U n is said to be O a . s . (v n ) (resp. o a . s . (v n )) if U n /v n is almost surely 
bounded (resp. if U n /v n converges towards with probability 1). 

The statistical model corresponds to a discretized version X = (X(i/n)) i=1 n of a 
locally self-similar Gaussian process defined by (1). 

One of the ideas of our method is to construct some estimators by using some prop- 
erties of the variance of the increments of X or the variance of the increments of order 2 
of X. While considering the increments of X is conventional since the associated se- 
quence is stationary, considering the increments of order 2 (or of a higher order) could 
be stranger. However, the main interest relies upon the fact that the observations of 
the latter resulting sequences are less correlated than those of the simple increments' 
sequence. All these vectors can actually be seen as special discrete filtering of the vector 
X. Let us now specify some general background on discrete filtering and its consequence 
on the correlation structure. The vector a is a filter of length £ + 1 and of order v > 1 
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with real components if 

e i 
y~] ct*a g = 0, for j = 0, . . . , v — 1 and ^ g^Og / 0. 

5=0 q=0 

For example, a = (1, —1) (resp. a = (1, —2, 1)) is a filter with order 1 (resp. 2). Let X a 
be the series obtained by filtering X with a, then: 

*°(£)=x>,*(^) fa^+i. 

Applying in turn the filter a = (1, — 1) and a = (1, — 2, 1) leads to the increments of X, 
respectively the increments of X of order 2. One may also consider other filters such as 
Daubechies wavelet filters, e.g. Daubechies (1992). 

The following assumption is needed by different results presented hereafter: 
Assumption Ai(k) : for i = 1, . . . , k 

v ®(t) = a 2 /3(i)|t| 2 ^- i + o (|i| 2H_i ) 

with (3(i) = 2H(2H - 1) . . . (2H -i + l) (where k > 1 is an integer). 

This assumption assures that the variance function v(-) is sufficiently smooth around 0. 
It allows us to assert that the correlation structure of a locally self-similar discretized 
and filtered Gaussian process can be compared to the one of the fractional Brownian 
motion. This is announced more precisely in the following Lemma. 

Lemma 1 (e.g. Kent and Wood (1997)) Let a and a' be two filters of length £ + 1 
and £' + 1, of order v and v' > 1. Then we have: 

= 7r'(i)(i + C°'(i)), (2) 

7n' a '(j) = 4r7 a ' a '(j), 7 a,a '(j) = ~\ £ a M« - J + ( 3 ) 



with 



and 



n 2H < w,> - w 2 

g,g'=0 



^Wk-«' + j| 2ff xr(^) 



°n U) 7 a ' a '(i) ' 
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Moreover, as \j\ — ► +00 

^ a '^)=o( W ^)- (5) 
Finally, under Assumption A\{u + i/), as n — > +cx) 

C°'(j)=o(l). (6) 

Remark 1 In i/ie case o/ i/ie fractional Brownian motion the sequence S n is equal to 0, 
whereas it converges towards for more general locally self-similar Gaussian processes, 
such as the Gaussian processes with stationary increments and with variance function 
v(t) = 1 — exp(— \t\ 2H ) or v(t) = log(l + \t\ 2H ) for which Assumption A\(k) is satisfied 
(for every k > 1). 

Remark 2 The stationarity of the increments and the local self- similarity required on 
the process X{-) are important, if the process is observed at times i/n for i = 1, . . . ,n. 
The crucial result of Lemma 1 is that the variance function of the filtered series behaves 
asymptotically as 7°(0). It seems to be difficult to relax the constraint of stationarity. 
Consider for example the Riemann-Liouville's process, e.g. Aids et al. (1999). This 
process is a Gaussian process which is H-self similar Gaussian but with increments sat- 
isfying only some kind of local stationarity. Following the computations of Lim (2001), 
the variance of the increments ' series of the Riemann-Liouville 's process is equal to 

2 N 




with I = Jq ((1 + u)^ -1 / 2 — u^ 1 / 2 ) 3 du + Jq U u 2H ~ 1 du. This integral cannot be asymp- 
totically independent of time. Note that this could be the case if the process is observed 
at irregular subdivisions. This question has not been investigated. 

Define Y a as the normalized vector X a with variance 1. The covariance between 
Y a (i/n) and Y a ' (i + j/n) is denoted by pn' a (j). Under Assumption A\{y + u'), the 
following equivalence holds as n — > +oo 

n a > a (j) ru p a > a (j) = ' (7) 

When a = a', we set, for the sake of simplicity 7™(-) = Jn' a (-), $%{■) = o~n' a (-), p%( ) = 
pT(-), 7 a ' a (-) = 7 a (0 and p a {-) = p^(-). 
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3 New estimators of H 

3.1 Estimators based on a convex combination of sample quantiles 

Let (p, c) = (pk,Ck)k=i,...,K G ((0,1) x for an integer 1 < K < +oo. Define the 

following statistics based on a convex combination of sample quantiles: 

K 

a P ,c;X a ) = Y J C k Z(Pk;X a ), (8) 
k=i 

where Ck, k = 1, . . . , K are positive real numbers such that X^fcLi c k = 1- For example, 
this corresponds to the sample median when K = l,p = l/2,c = 1 , to a mean of quar- 
tiles when K = 2,p = (1/4, 3/4), c = (1/2,1/2) . Consider the following computation: 
from Lemma 1, we have, asm +oo 

2 

Z(p,c;X a ) ~ 7 °(0)e (p,c;Y a ) . 

Remark 3 It may 6e expected that £ (p, c; Y a ) converges towards a constant as n — > 
+oo. /n ifee//, t/iis result is not interesting, since two parameters remain unknown: a 2 
and H and thus, it is impossible to derive an estimator of H. 

Remark 3 suggests that we have to use at least two filters. Among all available filters, 
let us consider the sequence (a m ) m >i defined by 



a? 



aj if i = jm 

for i = 0, . . . , ml , 

otherwise 



which is none other than the filter a dilated m times. For example, if the filter a = a 1 
corresponds to the filter (1,-2,1), then a 2 = (1,0,-2,0,1), a 3 = (1,0,0,-2,0,0,1), 
. . . As noted by Kent and Wood (1997) or Istas and Lang (1997), the filter a m , of length 
ml + 1, is of order v and has the following interesting property : 

7 am (0) =m 2H 1 a (0). (9) 

From Lemma 1, this simply means that E (X a? "(i/n) 2 ) = m 2H F, (X a (i/n) 2 ) , exhibit- 
ing some kind of self-similarity property of the filtered coefficients. As specified in the 
introduction, the same property can be pointed out in the context of wavelet decompo- 
sition. 
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Our methods, that exploit the nice property (9), are based on a convex combination 
of sample quantiles £ (p, c; g(X am )) for two positive functions g(-): g(-) = \ ■ \ a for 
a > and g(-) = log | • |. For such functions g(-) we manage, by using some property 
established in Lemma 1, to define some very simple estimators of the Hurst exponent 
through a simple linear regression. Other choices of the function g(-) have not been 
investigated in this paper. At this stage, let us specify that our methods extend the one 
proposed by Stoev et al. (2006); indeed they only consider the statistic £ (p, c; g(X am )} 
for p = 1/2, c = 1, g(-) = (-) 2 , that is the sample median of the squared coefficients. 
From (3) and (9), we have 

f (p, c; \X« m D = E ((X« m (l/n) f) a/2 f (p, c; |V» m |«) 



™ aH ^n^T' 2 (1 + C(0)) a/2 ^(P, C ;I^D, (10) 



n 

and 

1 



£(p,c;log|X a |)=-logE(X a (1/n)) 2 +£(p,c; log |) 



= tflog(m)+log( ^F7 a (0) 



a 2 



+ i log (1 + Sf (0)) + f (p, c; log |F- m |) . (11) 

Denote by = n~ 2 ' f/ (J 2 7 a (0). Equations (10) and (11) can be rewritten as 

log£ (p,c; |X am r) = afriog(m) +log (^ /2 C|y|- (p,c)) + e° , (12) 
^(p, C ;log|X am |) = Flog(m)+log(^) + 6og |Y| (p, c) + ^ g (13) 

with the random variables and £m E respectively defined by 



log (ikl^^) + " lo g(l + 5f f (0)), (14) 



and 

# = Up, c; log | r™ |) - £ log |y| (P, c) + 1 log (l + 5f (0)) (15) 

where, for some random variable Z, £z (p, c) = X)fcLi c& iz{pk)- We decide to rewrite 
Equations (10) and (11) as (12) and (13), since we expect that e m and e m g converge 
(almost surely) towards as n — > +oo. 

From Remark 3, two estimators of H can be defined through a linear regression 
of (logf (p,c; | I flm | a )) m=1> ... iM and (£ (p, c; log |) ) m=1> M on (logm) m=1) M 
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for some M > 2. These estimators are denoted by H a and H log . By denoting A the 
vector of length M with components A m = \ogm — jj J2m=i l°g( m )> m = I,. . . , M, we 
have explicitly from (12) and (13) and the definition of least squares estimates (see e.g. 
Antoniadis et al. (1992)): 

^=^M^i x °™i">L, „■ < 16 > 

^^(^•^i^L, m ' < 17 > 

/ j \ 1/2 

where \ \z\ \ for some vector z of length d denotes the norm defined by ( Yli=i z i ) 
We can point out that H a and H log are independent of the scaling coefficient a 2 . 



3.2 Estimators based on trimmed means 



Let < (3\ < (3<i < 1 and (3 = /?2), denote by g(X a ) the (3— trimmed mean of the 
vector g(X a ) given by 

1 ^ L M1J [n/3i]+l 

where {g(X a )) {l)n < (g(X a )) (2)n < ... < {g{X a )) {2)n are the order statistics of 
, ■ ■ ■ , (0(^ a ))„- It is well-known that {g{X a ))^ n = £(±;g(X a )). Hence, 
by following the ideas of the previous section, one may obtain 

log (jX° ro |« (/3) ) = aH log(m) + log (k^ /2 ^" (/3) ) + C* m , (18) 
log|X« m | (/3) = log(m) + log(K H ) + bg^ (/3) + e^'* m (19) 



with 



and 



e^' tm = log \Y" m \ (0) - log |y| (/3) + \ log (1 + , (21) 



where for some random variable Z, Z^ is referring to 



Z(P) = i f ^ ^(P)*- ( 22 ^ 

1 - P2 - Pi Jfi x 
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As in the previous section, two estimators of H, denoted by H a > tm and H log > tm , is derived 
through a log-linear regression 

H a ' tm = (\Xa m \<* {f3) ) . (23) 

a||A|| 2 V /m=l,...,M 

£log,ta = (l og \ X a m \ {/3) ) . (24) 

||A|| 2 V /m=l,...,M 

Remark 4 The estimator referred to the "estimator based on the quadratic variations" 
in the simulation study and studied with the same formalism by Coeurjolly ( 2001 ) cor- 
responds to the estimator H a,tm with a = 2, f3\ = fa = 0. 



4 Main results 

To simplify the presentation of different results, consider the two following assumptions 
on different parameters involved in the estimation procedures 

Assumption A2(p,c) : a is a filter of order v > 1, a is a positive real number, p 
(resp. c) is a vector of length K (for some 1 < K < +oo) such that < pk < 1 (resp. 
Cfc > and Ylk=i c k = 1)> M is an integer > 2. 

Assumption As((3) : a is a filter of order v > 1, a is a positive real number, 
f3 = (Pi,f3 2 ) is such that < (3\ < fa < 1, M is an integer > 2. 

Since A T (log(m)) m=1> M = \ \A\\ 2 and A T 1 = (where 1 = (l) m =i,..,M), we have 

A T ~, A T , 

H a_ H = e « and # lo s — H = e l ° s , (25) 

a||A||^ ll-^ll 

and 

ff Q ' im - # = - - e Q - tm and H log ' tm - H = — — £ ios ' tm , (26) 
a||A|| 2 \\A\\ 2 

wnere e — \£ m ) m =\,...,M' £ \ £ ™ )m=i,...,M> e ~ ^ Ern >m=i,...,M anct £ ~ 
{em' tm ) m= i m- Hence, in order to study the convergence of different estimators, it is 
sufficient to obtain some convergence results of sample quantiles £(p, g(Y a )) for some 
function g(-) and some filter a. Therefore, we first establish a Bahadur representation 
of sample quantiles for some non-linear function of Gaussian sequences with correlation 
function decreasing as k~ a , for some a > 0. In fact, the existing litterature on nonlinear 
function of Gaussian sequences (e.g. Taqqu (1977)) allows us to slighlty extend this 
framework by considering correlation function decreasing as k~ a L(k), for some slowly 
varying function L(-). 
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4.1 Bahadur representation of sample quantiles 

Let us recall some important definitions on Hermite polynomials. The j-ih Hermitc 
polynomial (for j > 0) is defined for t £ R by 

The Hermite polynomials form an orthogonal system for the Gaussian measure. More 
precisely, we have E (Hj(Y)Hk(Y)) = j\ 5j t k- For a measurable function g(-) defined on 
R for which E(g(Y) 2 ) < +oo, the following expansion holds 

9®=Y,% H &) with c i = ' E is(X)H j {Y)), 

where the integer r defined by r = inf {j > 0, Cj ^ 0}, is called the Hermite rank of the 
function g. Note that this integer plays an important role. For example, it is related to 
the correlation of g(Y\) and 9(^2) (for Y\ and Yi two standard gaussian variables with 
correlation p) since E(g(Y 1 )g(Y 2 )) = J2 k >r ^E~P k ^ P T \\9\\v^y 

In order to obtain a Bahadur representation (see e.g. Serfling (1980)), we have to 
ensure that F' g , Y \(£(p)) > and F' g ', Y J-) exists and is bounded in a neighborhood of 
£(p). This is achieved if the function </(•) satisfies the following assumption (see e.g. 
Dacunha-Castelle and Duflo (1982), p. 33). 

Assumption A4(£(p)) : there exist Ui, i = 1, . . . , L, disjoint open sets such that 
Ui contains a unique solution to the equation g(t) = i g (Y){p)i such that F^ Y ^{i{p)) > 
and such that g is a C 2 — diffeomorphism on L)^ =1 Ui. 
Note that under this assumption 

F g{Y) ^ g{Y) (p)) = f g (Y)(Zg(Y)(p)) = 2^ g>{g -l { ^ p)) y 



1= 



where gi(-) is the restriction of g(-) on Ui. Now, define, for some real u, the function 
MO by: 

h u (t) = l {gmu }(t)-F g{Y) (u). (28) 

We denote by r(n) the Hermite rank of h u (-). For the sake of simplicity, we set t p = 
T (£g(Y){p))- For some function g(-) satisfying Assumption A 4 (£(p)), we denote by 

r p = inf r( 7 ), (29) 

7€Uf =l9 (i7 i ) 

that is the minimal Hermite rank of h u (-) for u in a neighborhood of £ g (Y)(p)- 
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Theorem 2 Let {Y(i)}^~J^ be a stationary (centered) gaussian process with variance 1, 
and correlation function p(-) such that, as i — > +oo 

\p(i)\ ~L(i)r Q , (30) 

/or some a > and some slowly varying function at infinity L(s),s > 0. T/ten, under 
Assumption A4(£(p)), we have almost surely, as n — ► +oo 

Z(P;9(Y)) -Zg(Y)(P) = 7 77 7-77 + Oa.s. (r„(a,T p )), (31) 

Jg(Y)KSg(Y)\P)) 

the sequence (r n (a,T p )) n>1 being defined by 



r n {a,T p ) = < 



(32) 



' n~ 3 / 4 log(n) 3 / 4 if ar p > 1, 

n - 3 / 4 log(n) 3 / 4 L Tp (n) 3 / 4 i/ ar p = 1, 

n -l/2-ar p /4 log ( n )r p /4+l/2 L ( n )r p /4 ^ 3/3 < ^ < ^ 

n - a ~p log(n)~f L(n)~p i/ < af„ < 2/3, 

where for some r > 1, L r (n) = X)|i|< n l/°(*)| T - 

Note that if L(-) is an increasing function, L T (n) = (log(n)L(n) T ). 

Remark 5 Without giving any details here, let us say that the behaviour of the se- 
quence r n (-,-) is related to the characteristic (short-range or long-range dependence) of 
the process {h u (Y (i))}f =1 for u in a neighborhood of £ g (y)(p). In the case af p > 1, 
corresponding to short-range dependent processes, the result is similar to the one proved 
by Bahadur, see e.g. Serfling (1980), in the i.i.d. case. For short-range dependent 
linear processes, using a law of iterated logarithm's type result Wu (2005) obtained a 
sharper bound, that is n~ 3 / 4 loglog(n) 3 / 4 . This bound is obtained under the assump- 
tion that F'(-) and F"(-) exist and are uniformly bounded. For long-range dependent 
processes (aT p < 1), we can observe that the rate of convergence is always lower than 
n~ 3 / 4 log(n) 3 / 4 and that the dominant term n~ 3 / 4 is obtained when ar p — > 1. 

We now propose a uniform Bahadur type representation of sample quantiles. Such 
a representation has an application in the study of trimmed-mean. For < po < p\ < 1 
consider the following assumption which extends A4(£(p)) 
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Assumption A^(po,pi) : there exists Ui, i = 1, . . . , L, disjoint open sets such that 
Ui contains a solution to the equation g(t) = £ g (Y)(p) f° r a h Po < P < Pi, such that 
F' g pr)(ti(p)) > for all po < p < p\ and such that g is a C 2 — diffeomorphism on \jf =1 Ui. 
Under the previous assumption, define 



'Pa, Pi 



inf t(7). 

ieuf =1 g(Ui) 



(33) 



Theorem 3 Under the conditions of Theorem 2 and Assumption As(po 5 Pi); we have 
almost surely, as n — > +oo 



sup 

po<p<pi 



i(p;g{Y))-i g{Y) (p) 



fg(Y)(tg(Y)(p)) 



O a .s. (r n {a ))• (34) 



Remark 6 To obtain convergence results of estimators of H , some results are needed 
concerning sample quantiles of the form £ (p; g(Y am )) , with g(-) = \ ■ |. Lemma 14 
asserts that the Hermite rank t p of the function ^(-) with g(-) = \ ■ \, is equal 
to 2 for all < p < 1. Moreover, for all < p < 1 and for all < po < p\ < 1, 
Assumptions A4(£(p)) and A^{po,pi) are satisfied, and we have T p = T P(hPl = 2. 
Since from Lemma 1, the correlation function ofY am satisfies (30) with a = 2v — 2H 
and L(-) = 1, by applying Theorem 2, the sequence r n (-, •) is then given by 



r n {2u - 2H, 2) = n~ 3 / 4 log(ra) 3 / 4 , if v > 2 



(35) 



and for v = 1 



r n (2-2H,2) = { 



n -3/4 i og ( n )3/4 ifO<H <3/4, 

n -3/4 log ( n )3/2 { f H = 3 /^ 

n -i/2-(i-f0 log(n) t/ 3/4 < H < 5/6, 

„-2(2-2f0 log(n) 2 »/ 5/6 < iJ < 1. 



(36) 



4.2 Convergence results of estimators of H 

In order to specify convergence results, we make the following assumption concerning 
the remainder term of the variance function v(-). 

Assumption Ae(n) : there exists rj > such that v(t) = a 2 \t\ 2H (1 + O (|i|' ? )) , as 
\t\ -► 0. 

The first result concentrates itself on estimators H a and H log based on a convex 
combination of sample quantiles. 
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Theorem 4 Under Assumptions A\(2u), A^ijp, c) and Aq(tj), 

(i) we have almost surely, as n — > +00 

O (n"") + O a ., K 1 / 2 log(n)) if v > H + \, 
H a -H={ 0{n-v) + O a . s .{n- l l 2 \og{nf/*) ifv = \,H = \, (37) 
(n"") + Oa. s . (n" 2 * 1 -^ log(n)) i/ 1/ = 1, § < # < 1. 

yl similar result holds for H log . 

(ii) the mean squared errors (MSE) of H a satisfies 

MSE (h° -H} =0 (v n (2u - 2H)) + O (r n (2u - 2H, 2) 2 ) + O (n~ 2r ') . (38) 

The sequence r n (2v — 2H,2) is given by (35) and (36) and the sequence v n (-) is defined 
by 

' n- 1 ifv>H+\, 
n- 1 log(n) ifv = l,H = l (39) 
„-4(i-H) if u = l,l<H<l. 

Again, the same result holds for MSE 

(m) i/ i/ie /iZter a is such that u > H + 1/4, and if rj > 1/2, then we have the 
following convergence in distribution, as n — > +00 

^n~(H a -H^j ^M(0,a 2 a ) and ^ (# log - tf) — ► Af(0, a 2 ), (40) 

where u 2 is defined for a > by 



v n {2v - 2H) = < 



r 



The vector B is defined by B = , and the real numbers qk and ir^ are defined by 



-i( l +Pk\ „„ A „ a _ (Qk 



Qk = *" x ( t-hP- ) and tf = -™L—. (42) 



Finally, the matrix R(i,j), defined for i e Z and j > \, is a M x M matrix whose 
(mi, 7712) enin/ is 

(^,i)) mi ^=p 0roi ' am2 W 2 ^ (43) 
where p ami > am2 (.) £5 ^ e correlation function defined by (7). 
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Remark 7 The expression of the variance a 2 a given by (41) could appear to be very 
complicated. However, given some vectors p and c and some integer M, it does not take 
unreasonnable effort to compute it for each value of H by truncating the two series. This 
issue is investigated in Section 5 to compare the different parameters. 

Remark 8 Let us discuss the result (38). The first term, O (v n ), is due to the variance 
of the sample cumulative distribution function. The second term, O (r^) is due to the 
departure of £ (p) — from F — p. We leave the reader to check that 

[ 0{v n (2v-2H)) %fv>H+\, 
0(r n (2u-2H,2f)+0(v n (2u-2H)) = { ) n{ " 

I 0(r n {2v-2H,2f) ifv<H+\. 

Finally, the third one, O (n _2r? ) is a bias term due to the misspecification of the variance 
function v(-) around 0. 

Remark 9 If K = 1, we have, for every a > 0, 

iez j>i H v •"' 

Assume ^6(77) with i] > 1/2 which allows to neglict the bias term with respect to the 
variance one. The result (40) is proved by using some general central limit theorem 
obtained in this dependence context by Arcones (1994), which is available as soon as 
p a (-) 2 is summable. Therefore, if only Ai(2) is assumed, the filter a cannot exceed 1 
(and then correspond to a = (1,-1)) and, due to (5), the result (40) is valid only for 
< H < 3/4. As a practical point of view, one observes that for such a filter and large 
values of H , the estimators have very big variance. Note that if Ai(2v) can be assumed 
for v > 1, then the asymptotic normality is valid for all the values of H. 

The next result asserts the link between H log and H a . 

Corollary 5 Let (a n ) n >i be a sequence such that a n — > 0, as n — > +00. Then, under 
conditions of Theorem 4 (H), the following convergence in distribution holds, as n — > +00 

V^(h^-h) -^M(0,a 2 ). (44) 

The following theorem presents the analog results obtained for the estimators H a < tm 
and H log,tm based on trimmed-means. 
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Theorem 6 Under Assumptions A\(2v), A^{fH) and Aq(t)), properties (i) and (ii) of 
Theorem 4 hold for the estimator H a ' tm and H log,tm with the same rates of convergences. 

(Hi) if the filter a is such that v > H + 1/4 and if n > 1/2, then, under the notations 
of Theorem 4, we have the following convergence in distribution, as n — > +oo 

(H a > tm -H)-^ M(0, a% tm ) and ^ (h 1 ° sMi - h) M(0, a 2 tm ), (45) 

where o~ a £ m is defined for ol > by 

withq=<$>- 1 (±±£). 

5 Numerical computation and simulations 

5.1 Asymptotic constants a 2 a and (T 2 a tm 

In order to compare the different estimators, we intend to compute the asymptotic con- 
stants <7 2 and tm defined by (41) and (46) for various set of parameters (a,p, c, (3, M). 
For this work, both series defining a\ and o" ft [ m are truncated (|z| 200, j ^ 150). Fig- 
ure 2 illustrates a part of this work. We can propose the following general remarks: 
• Among all filters tested, the best one seems to be 



a 



incl if0<#<3/4, 
dbA otherwise. 



where incl and db4 respectively denote the filter (1,-1) and the Daubechies wavelet 
filter with two zero moments explicitly given by 

dbA = (0.4829629,-0.8365763,0.22414386,0.12940952) . 

• Choice of M: increasing M seems to reduce the asymptotic constant u 2 . Obviously, 
a too large M increases the bias since £ (p, c; g(X aM or g(X aM ) are estimated 
with N — Mi observations. We recommend setting it to the value 5. 

• We did not manage (theoretically and numerically since series defining (41) and (46) 
are truncated) to determine the optimal value of a. However, for examples considered, 
it should be near the value 2. 
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• Again, this is quite difficult to know theoretically and numerically which choice of p 
is optimal. What we observed is that, for fixed parameters a, M and a, the asymptotic 
constants are very close to each other. 

• Choice of p in the case of a single quantile (see Figure 2): the optimal p seems to 
be near the value 90%. However, p = 1/2, corresponding to the estimator based on the 
median, leads to good results. 

• Choice of /3i = /?2 = /? for the estimators based on trimmed- means (see Figure 2): 
obviously the constant grows with (3 but we can point out that estimators based on 
10%— trimmed-means are very competitive with the ones obtained by quadratic varia- 
tions (ft = 0). 

5.2 Simulation 

A short simulation study is proposed in Table 1 and Figure 1 for n = 1000 and H = 0.8. 
We consider two locally self-similar Gaussian processes whose variance functions are in 
turn v(t) = \t\ 2H (fractional Brownian motion) and v(t) = 1 — exp(—\t\ 2H ). To generate 
sample paths discretized over a grid [0, 1], we use the method of circulant matrix (see 
Wood and Chan (1994)), which is particularly fast, even for large sample sizes. Various 
versions of estimators are considered and compared with classical ones, that is the one 
based on quadratic variations, Coeurjolly (2001), and the Whittle estimator, Beran 
(1994). In order to illustrate the robustness of our estimators, we also applied them to 
contaminated version of sample path processes. We obtain a new sample path discretized 
at times ijn and denoted by X c (i/n) for i = 1, . . . , n through the following model 

X c (i/n) =X(i/n) + U(i)V(i), (47) 

where U(i), i = 1, ...,n are Bernoulli independent variables #(0.005), and V(i), i = 
1, . . . ,n are independent centered Gaussian variables with variance such that the 

signal noise ratio at time i/n is equal to 20 dB. As a general conclusion of Table 1, 
one can say that all versions of our estimators are very competitive with classical ones 
when the processes are observed without contamination and they seem to be particularly 
robust to additive outliers. Both bias and variance are approximately unchanged. This 
is clearly not the case for classical estimators. Indeed, concerning quadratic variations' 
method, the estimation procedure is based on the estimation of E((X am (1/n)) 2 ) by 
sample mean of order 2 of (X ™) 2 , Coeurjolly (2001)), that is particularly sensitive to 
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additive outliers. Bad results of Whittle estimator can be explained by the fact that 
maximum likelihood methods are also non-robust methods. 



6 Proofs 

We denote by || • Ui^^) (resp. || • \\g q ) the norm defined by ||/i||l2(^ = E(h(Y) 2 ) 1 / 2 
for some measurable function h(-) (resp. (Yliez l^l 9 ) 1 ^ 2 f° r some sequence (ui)i & z). In 
order to simplify the presentation of proofs, we use the notations F(-), £(•), /(•), F (•) and 
£(■) instead of F g(Y) (-), £ ff (y)(-), fg(Y){'), F g (Y) (sf/OO) and (-,g(Y)) respectively. 
For some real x, [x] denotes the integer part of x. Finally, A denotes a generic positive 
constant. 



6.1 Sketch of the proof of Theorem 2 

We give here a brief explanation of the strategy to prove Theorem 2. This proof follows 
exactly the one proposed by Serfling (1980) in the i.i.d. case. One starts by writing 

p — F (£(p)) 



urn) 



(Z(p)-Z(pj)=A(p)+B(p) + C(p), 



with 



p-f (Up)) 

m = -^rnr (48) 

F ((<!>)) - F («p)) - (F({ (p)) - F(((p))) 

B <^^ w (49) 

cw=3i w M -(* , - £W )' (50) 

From the definition of sample quantile, we have almost surely, see e.g. Serfling (1980), 
A(p) = O a , s , (n^ 1 ). Now, in order to control the term C(p), Taylor's Theorem is used 
and a control of £ (p) — £(p) is needed. The latter one is done by Lemma 10 which 
exhibits the sequence e n (a,T p ) such that £ (p) — £(p) = O a , s , (e n (a, r p )). Then, in order 
to control B(p) it is sufficient to control the random variable 

S n (t( P ),e n (a,T p ))= sup |A(f(p) + a;)-A(£(p))|, 

\x\<S n (a,T p ) 

with A(-) = F (•) — F(-). This result is detailed in Lemma 11. In order to specify the rate 
explicited by Theorem 2, we present and prove Lemmas 10 and 11. Some preliminary 
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results, given by Lemma 7, Corollary 8 and Lemma 9, are needed. Among other things, 
Lemma 7 and Corollary 8 propose some inequalities for controlling the sample mean of 
non- linear function of Gaussian sequences with correlation function satisfying (30) . 

6.2 Auxiliary Lemmas for the proof of Theorem 2 

Lemma 7 Let {Y(i)}fJ^ a gaussian stationary process with variance 1 and correlation 
function p(-) such that, as i — > +oo, \p(i)\ ~ L(i)i~ a , for some a > and some slowly 
varying function at infinity L(-). Let h(-) G L 2 (d(f>) and denote by r its Hermite rank. 
Define 

1 n 
n r— 1 

Then, for all 7 > 0, there exists a positive constant (a, t), such that 

P(|F n | >^y n ) = 0{pT l ), (51) 



y n = y n (a,T) = < 



n 1 / 2 log(n) 1 / 2 if ar > 1, 

n-VaiogfajVa-Lr^Va */ ar = 1, (52) 

n -ar/2 l og ( n )T/2 L ( n )r/2 -J < «T < 1. 

where L T (n) = Yl\i\< n \p(^)\ T • ^ n the case ar = 1, we assume that for all j > r, the 
limit, \\m. n ^ +00 L T (nY l Y J \i\<n\p( i )\ j exists. 

Proof. Let (y n )n>i be the sequence defined by (52). The proof is splitted into three 
parts according to the value of ar. 

Case (xt < 1 : From Chebyshev's inequality, we have for all q > 1 



{\Y n \>K iyn )<^^^({Y n ) 2ll \ 

^7 Vn V / 



fc-y y n 

From Theorem 1 of Breuer and Major (1983) and in particular Equation (2.6), we have, 
asm +oo 

E (< 7 »> 2 ")~w^ w,th " 2 =EE¥"« J . (»> 

where Cj denotes the j-th Hermite coefficient of h(-). Note that a 2 < 1 1^1 1^2(^1 |p|||t- 
Thus, for n large enough, we have 

p(|F„| > K 72/n ) < (iWlWMlK 2 )*- (54) 
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From Stirling's formula, we have as q — > +00 

M „ V2q«(2e- l y. (55) 
From (52) by choosing q = [log(n)], (54) becomes 

P (\Y n \ > ^y n ) < A (2 e - 1 ||/ l ||| 2W) ||/ 3 ||i^ 2 ) 1 ° S(n) = O (n-~<) , 
if 4 >2||/i|| 2 L2(#) || /0 ||^exp( 7 -l). 

Case olt = 1 : Using the proof of Theorem 1' of Breuer and Major (1983), we can prove 
that for all q > 1 

E ((nV% (n) -V2 Fn) 2 9 ) < A ^_ E ( (n i/2 Lr(n) -i/2 Fn)2 y 



iy 3 



<x p}(y(^l_ lim Lr{n) -iy 1 ( 

H \j>r J \i\<n 

<^M% m - (56) 
Then from Chebyshev's inequality, we have for all q > 1 

r(|F„| <^M(| W | W - 2 )». 

From (52) by choosing q = [log(n)], we obtain 

IP (\Y n \ > K jyn ) < A (2e- 1 ||/ i ||i 2(#) K- 2 ) 10 ^ = O (n-T) , 
if 4 > 2\\h\\\ 2m x exp( 7 - 1). 

Case cur < 1 : Denote by k a the lowest integer satisfying k a a > 1, that is k a = [l/a]+l, 
and for j > r denote by Zj the following random variable 

1 n 

Denote by «i j7 and K2 i7 two positive constants such that k 7 = max(/«i )7 , K2, 7 )- From the 
triangle inequality, 

P (\Y n \ > K lVr ) < P(jF n -J2Zj\> ki, 7 2/„) + J2¥(\Zj\> n 2nVn ) (57) 
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Since 

k(y. 1 



j=r i=l j>k a i=l 

where h'(-) is a function with Hermite rank k a . Applying Lemma 7 in the case ar > 1, 
it follows that, for all 7 > 0, there exists a constant K\ a such that, for n large enough 

k — 1 

P(|F„ - |] Zj| > Ki )7 j/ n ) = O (n"T) . (58) 
Now, let t < j < k a and 5 > 1, from Theorem 3 of Taqqu (1977), we have 

p(|^-|>«2 )7 y„) <^-^-A%) V 2 «e( ^ ^ ( y(n))...^(y(^)) 

K 2j7 2/n \3-J \h,...,i 2q J 



L(n\"i fcj _^ 2q 



where is a constant such that fi2 q < (iz^j) E (Hj(Y) 2q ) . It is also proved in Taqqu 
(1977) (p. 228), that E {H^Y) 2 *) ~ (2jq)\/(Z>*(jq)\), as q -> +00. Thus, from Stirling's 
formula, we obtain as q — ► +00 



(l Z il^«)<A^ r lo g („)-V^ T f-^j (^j «Sj • 
By choosing g = [log(n)], we finally obtain, asm +00 

fc °~ 1 / /, \ 2/9 T \ T \log(n) 

E P (l^l^^n)<A(^ r -^(^) (-J k^J =0(n^), (60) 



if /«2 i7 > 1 _ 2 Qr (fy) (2r) T exp(7 — r). From (57), we get the result by combining (58) 
and (60). ■ 



Corollary 8 Under conditions of Lemma 7, for all a > 0, j > 1 and 7 > 0, there exists 
q = (7(7) > 1 and C7 > swc/i £/iat 

E ^E^-( y w)} j <c 7 ^ 7 - (si) 
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Proof. (53), (56) and (59) imply that there exists A = A(g) > such that for all q > 1, 
we have 



E ^E^( y «)} j<%)n^ = A(g)x< 



if aj > 1 



L Tp (n)n q if aj = 1 
L(n) ajq n~ a i q if aj < 1 



= O (n~ 7 ) . (62) 
Indeed, it is sufficient to choose q such that, q > 7 if aj > 1 and g > 7/07 if aj < 1. 



Lemma 9 Lei < p < 1, denote by g(-) a function satisfying Assumption A4(£(p)) 
and by (x n ) n >i a sequence with real components, such that x n — ► 0, as n — > +00. T/ien, 
/or a// i > 1, i/iere existe a positive constant dj = dj(£(p)) < +00 such that, for n large 
enough 

\cj{€(p) + x n ) - Cj(£(p))\ < dj \x n \. (63) 

Proof. Let j > 1, under Assumption A4(£(p)), for n large enough, + x n £ 
Ui =1 g(Ui). Thus, for n large enough, 

cMp) + *n) ~ cMp)) = [ (h( p)+Xn (t) - h m (t)) Hj(t)cf>(t)dt 

JR 
L 

= ( 1 9,W<5(p)+^n - 1 » i (t)<€(p)) H j{t)<t>{t)dt 

i=i ,/c/ » 



( 



Eti - ( W,n) - 0(mi,„)) if J = 1, 

k Eil(-l) J '(^'- 1) (M,n)-^'- 1) (mi 1 „)) if i > 1, 

where ^(-) is the restriction of #(■) to C/j, and where mj jra (resp. Mj j7l ) is the minimum 
(resp. maximum) between g~ l {i{p) + x n ) and g^ 1 (^(p))- We leave the reader to check 
that there exists a positive constant dj, such that, for n large enough 



^-^1 = 1 

|cj(£(i>) +x n ) -Cj(£(p))\ <dj\x n \ x <( 

x Ej=i 

which is the desired result. ■ 



^(gt'Hu)) if J = 1,2 

^- 2 )(^- 1) (u)) if J > 2, 
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Lemma 10 Under conditions of Theorem 2, there exists a constant denoted by k £ = 
K £ (a,T p ), such that, we have almost surely, as n — > +oo 



t(p;g(Y))-t g ( Y )(p) 



< Sr. 



where e n = s n (a,r(^{p))) = K £ y n (a,T(£(p)), y„(-,-) being defined by (52). 
Proof. We have 

p (|f(p) - e(p)| > £«) = p (f(p) < e(p) - en) + p (f (?) > + *«) . 

Using Lemma 1.1.4 (m) of Serfling (1980), we have 

P (f(p) < £(p) " en) < P (£(p) - s n ) > p) . 
Under Assumption A 4 (£(p)), for n large enough 

p - F(H(p) - e n ) = /(£(p))e„ + o {e n ) > ^P±e n . 



Consequently, for n large enough and from (66) 



(f (p) < £(p) - e n ) < P (f (£(p) - e n ) - F(£(p) - e n ) > 



m( P )) 



Define r Pjn = r(£(p) — e n ), from Lemma 9, we have for n large enough 

£ (£(P) " £n) " *U(p) " £n) > 2 (F (£(p)) - F(C(p))) + 2 £ „ Y, Z n, 



where 



= < 



(64) 



(65) 



(66) 



(67) 



(68) 



{Tp ^ 3 — Tp,n\ if Tp,n > Tp-, 

1 n d 

ifr Piri = r p , and Z n>j = - ^jH^Y {i)). 

{Tp,n ^ j ' < 7p} if 7p,n < 7p- 



Now, define c e = ft e /(£(p))/4. Let 7 > 0, (61) implies that there exists g > 1 such that, 
for n large enough 



P (|2e„Z n | > < X)P(|Z nj - 



* E i E (ZZ) = O (n-) 



(69) 
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Let us fix 7 = 2. From (67), (68) and (69) and from Lemma 7 (applied to the function 
h,£(p\(-)), we obtain 

P (f(p) < e(p) - £n) < P (|£ (£(p)) - *U(p))I > c £ £n) + O (n- 2 ) = O (n- 2 ) , 

if c £ > K2 that is if k £ > 4//(£(p))«2. 

Let us now focus on the second right-hand term of (65). Following the sketch of this 
proof, we may also obtain, for n large enough 

P(f(p)>£(p) + £n) =0{n~ 2 ), 

if k £ > 4//(^(p))fv2- Thus, for n large enoug h P (|£(p) - f(p)| > e n ) = O (n- 2 ), which 
leads to the result thanks to Borel-Cantelli's Lemma. ■ 

The following Lemma is an analogous result obtained by Bahadur in the i.i.d. frame- 
work, see Lemma E p. 97 of Serfling (1980). 

Lemma 11 Under conditions of Theorem 2, denote by A(z) for z £ R the random 
variable, A(z) = F (z; g(Y)) — F g (y)(z). Then, we have almost surely, as n — > +oo 

Sn(£,g(Y){p),£n{a,T p )) = SUp |A(£ ff (y)(p) + x) - A(£ g( y)(p))\ = O a . s . (r n (a,T p )) , 

|x|<£„ 

(70) 

where e n = e n (a,T p ) is defined by (64) and r n (a,T p ) is defined by (32). 

Proof. Put s n = e n (a,T p ) and r n = r n (a,T p ). Denote by (/3„)„>i and (jib,n)n>i the 
following two sequences 

Pn = [n 3/4 e n ] and rj b>n = £(p)+e n -^-, 

Pn 

for b = —P n , . . . , P n - Using the monotonicity of F(-) and F (•), we have, 

S n (Z(p),£n)< max \M b>n \ + G n , (71) 

-Pn<0<Pn 

where M b>n = A(r/ M ) - A(£(p)) and G n = max„ /3n <K /3 „_ 1 (F(rj b+1;n ) - F{r) b , n )) . Under 
Assumption A^^p)), we have for n large enough 

G n < (Vb+l,n - Tj b ,n) X SU P /(£(p) + x) = O (V 3 / 4 ) . (72) 

\x\<6„ V ' 
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The proof is finished if one can prove that for all 7 > (in particular 7 = 2) and for all 
b, there exists k'^ such that 

F(\M b:n \> K l J r n )=0(n-^). (73) 
Indeed, since (3 n = O (n 1//2+<5 ) for all 6 > 0, if (73) is true, then we have 

P( max \M b>n \ > K 2 r n {a,T p )) < (2/3 n + 1) x max P (|M 6 , n | > 4r n ) 

-Pn<0<Pn -p n <b<P„ 

= O (n-^ 2+& ) . 
Thus, from Borel-Cantelli's Lemma, we have, almost surely 

max \M b>n \ = O a . s . (r n ) 

-Pn<0<Pn 

And so, from (71) and (72). 

S n (Z(p),e n ) = O a . s . (r n ) + O (n- 3 / 4 ) = O a . s . (r n ) , (74) 

which is the stated result. 

So, the rest of the proof is devoted to prove (73). For the sake of simplicity, denote 
by h' n (-) the function h Vbn (-) — h^(-). For n large enough, the Hermite rank of h' n {-) is 
at least equal to T p , that is defined by (29). In the sequel, we need the following bound 

fOT \Wn\\h {d4>) 

\K\\l*W) = E ( /l n( y ) 2 ) = ^(1 ~ UJ n ) With Wn = \F g{Y) ( Vb , n ) - F g(y) (Z(p))\. 

As previously, we have u n = 0{e n ) and so, there exists ( > 0, such that 

\K\\h W) < Ce n . (75) 
From now on, in order to simplify the proof, we use the following upper-bound 

£ n £ n (o:, Tp) < e n (a, Tp), 

and with a slight abuse, we still denote e n = e n (a,T p ). Note also, that from Lemma 9, 
the j-th Hermite coefficient, for some j >t p , is given by Cj(r) bin ) — Cj(£(p)). And there 
exists a positive constant dj = dj(£(p)) such that for n large enough 

\Cj(Vb,n) ~ Cj{£(p)\ < dj EjM- < dj £ n . (76) 

Pn 
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We now proceed like in the proof of Lemma 7. 

Case OLTp > 1: using Theorem 1 of Breuer and Major (1983) and (54), we can obtain 
for all q > 1 

P (lAf^l > <r B ) < A^M_J_| K | Sw) || p || ?p . (77) 
As q — > +oo, we get 



P (|M M | > ^r n ) < A-^-?« ^Ce" 1 



1 



From (32), (52) (with r = r p ) and by choosing g = [log(n)], we have 

/ i \ log(n) 

P (|M 6 , n | > <r n ) < A ^C^e-^lpH^^J = Q (n^) , 
if K ; 2 >2CK e ||p||^exp( 7 -l). 

Case cuTp = 1 from (56), we can obtain for all q > 1 

E (Mf„) < A<M ^1\ K \\« < A <« M 

L= (n) q e q n , 
< A TpV ; (2C6" 1 ) 9 g«. 

From (32), (52) (with r = r p ), by choosing q = [log(n)], we have 

P(|M 6 , n | > <r n ) < * E [Mil) 



,2 , \ lo §( n ) 



(78) 



if k^ 2 > 2(K £ d 2 Tp /T p \ exp( 7 - 1). 

Case orFp < 1: denote by (ri >n ) n >i and by (r2, n ) n >i the following two sequences 

r 1>n = n -V2-aT P /4 log( n )^/ 4+1 / 2 L(n)^/ 4 and r 2 ,„ = n~ a ^ log(n) T *L(n) T »\ (79) 

Note that max (ri >n , r2, n ) is equal to ri >n , when 2/3 < orf p < 1 and to r2, n , when 
< crf p < 2/3. So, in order to obtain (73) in the case < af p < 1, it is sufficient to 
prove that there exists k'^ such that, for n large enough 

P (|M 6 , n | > K ;max(r lin ,r 2in )) = O (n^) . 
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Denote by k a the integer [I /a] + 1 for which ak a > 1, and by Zj >n for f v < j < k a the 
random variable defined by 

n i=l 

From the triangle inequality, we have 



IP (|M M | > K ;max(n, n ,r 2in )) < P(|M M - ^ Z j>n \ > ^r lin )+ ^ P (|%„| > . 
Since, 



j=Tp j=T p 

(80) 



3=T P t=l j>k a i=l 

where /i"(-) is a function with Hermite rank k a , such that a/c a > 1, we have from (77) 

P(|M M - p 3 , n | > <n >fl ) < A ^\K\\% {d ^^ q \\A ( 81 ) 

for all g > 1. From (75), we obtain, as q — > +cx) 

P(|M M - £ Z j)n | > «;r 1>n ) < A-^-g* (2Ce- 1 ||p||^<- 2 ) 9 . 

i=r p n r l,n 

From (52) (with r = r p ), (79) and by choosing q = [log(ra)], we obtain 



P(|M M - £ > «X„) < A (2Ce- 1 || /9 ||^ K£K ; 2 ) l0g(n) = O (n" 7 ) , (82) 



if k!^ > = 2(| |p| k e exp(7 — 1). Now, concerning the last term of (80), from (59), 
we can prove, for all r p < j < k a 

» \ L(n)i q 1 fcAribn) - c i (£(p))\ 2q 
P ^> > <r 2 , n < A V ; — jKlb ' nJ J ^ yP,) fi 2q , 



where is a constant such that, as q — > +oo, 

u <A ( 2 V (2jg)! 
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From (76), we have, as q — > +00 

{Zj, n > K r 2 , n ) < A : — ^—q> q - : — ) d , /e 

From (32), (52) (with r = r p ) by choosing q = [log(n)], we have, asm +00 



n a / \ 1 — aj \ e 



Consequently, as n — > +00, we finally obtain 

fe a — 1 / / _\ t \ l°g( n ) 

EPfe>W<A ^ ffl 4- £ 2 ^ 2 =0(n^), (83) 

i=r p V / 

if k! 2 > k' 2)7 = 1 2 a - (^") T 4 Kgexp(7 — r). Let us choose such that k! 2 > 
max(/«' 1)7 , ^2,7)- Then, by combining (82) and (83), we deduce from (80) that, for every 
7 > 

P (|M 6 , n | > ^max(r lin ,r 2)n )) = O (n" 7 ) , 
and so, (73) is proved. ■ 

6.3 Proof of Theorem 2 

Proof. Let us detail the proof presented in Section 6.1. We have 
p-F{£{p)) 



(i{p)-ttp))=A(p)+B(p) + C{p) 



with A(p), B(p) and C(p) respectively defined by (48), (49) and (50). Under Assumption 
A-4(£(p)) , from Lemma 10 and Taylor's theorem we have almost surely, as n — > +00 

C{p) < sup F; (y) (£(p) + s) (f (p) - £(p)) 2 = Oa. s . (e„(a, r p ) 2 ) . 

\x\<s„(a,T p ) 

From the definition of sample quantile, we have almost surely, see e.g. Serfling (1980), 
A(p) = O a , s , (fi -1 )- Now, by combining Lemma 10 and Lemma 11, we have almost 
surely B(p) = O a , s , (r n (a,T p )). Thus, we finally obtain 

Zip) ~ m = P ~ f ^f ] + + O a .s. (r„(a,r p )) + O a . a . {e n (a,r p ) 2 ) , 

which leads to the result by noticing that e n (a>, r p ) 2 = O (r n (a,r p )). ■ 
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6.4 Auxiliary Lemmas for the proof of Theorem 4 

Let < p < pi < 1. 

Lemma 12 Under conditions of Theorem 3, there exists a constant denoted by 6 = 
9(a,T POjPl ) such that, we have almost surely, as n — ► +oo 



T = sup 

Po<P<Pi 



), (84) 
w/iere e n = £ n (a, t POjP1 ) = 6y n (a,T POtPl ) and y n is given by (50). 



Proof. Define p,-^ = Po + t^st (Pi ~Po) for j = 0, . . . , [n 3/2 ], and let p G [p ,Pi]. Using 
the monotonicity of £(•) and £(•), there exists some j such that p G \pj,n,Pj+i,n] an d such 
that 

£ (p) - £(P) < £ (?) - f (Pi+l,n)) + f(Pi+l,n)) - £(p) 

< £ (Pj+l,n) - £(Pj+l,n)) + £(Pj+l,n)) - £(p,>)) + £(P?>)) ~ £(p) 

< £ (Pj+l,n) - £(Pj+l,n)) + £(Pj+l,n)) - £(Pj>))- 

This leads to 



T < max 

i=0,..,[n 3 / 2 ] 



£(Pj,n) -£(Pj,n)) + max |£(p,+l, n ) ~ £(P?>))I • (85) 
j^O,...,^/ 2 ]-! 



Under Assumption A 5 (p 5Pi)) it comes 

max |£(p j+ i,„) - £(p,>))| = O (V 3 / 2 ) . 

Now, following the proof of Lemma 10, one can prove that there exists some constant 
6(a, t POiP1 ) such that for all j = 0, . . . , [n 3 / 2 ], 

l£ (Pj,n) - £(Pj,n)| > 6>y n (a,T POiPl )) = O (n~ 3 ) . 
Therefore, as n — > +00, 

max |£ (pj,„) - £(pj,„)| >£„) < ([n 3/2 ] + 1) max P (|£ (Pj,n) - £(Pj,n)| > 

=0,...,[n J /^] / jM),...,^/ 2 ] ^ 

= O (n- 3 / 2 ) . 

which, combined with (85), (86) and Borel-Cantelli's Lemma, leads to the result. ■ 

The following result is an extension of Lemma 11 and Theorem 4.2 obtained by Sen 
(1971). 
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Lemma 13 Under Assumptions of Theorem 3 and following Lemma 11, we have almost 
surely, as n — > +oo 

S*= sup \A(x)-A(y)\ = O a . s .(r n (a,T P0 , Pl )) (87) 

x,y e [£(po),£(pi)1 

k — 2/1 < e n (a,T P0 . pi ) 
where t P0)Pi is defined by (33). 

Proof. Set e n = e n (a,T POtPl ) and r n = r n (a,T POtPl ). Define £ j>n = £(p ) + ^(£(pi) ~ 
£(po)) for j = 0, . . . ,p„ withp„ = [e" 1 ] , and let x, y G [£(po),£(pi)] such that < e„. 

Two cases may occur 

• If there exists some j such that x,y £ [Cj,n,^j+i,n] then 

|A(x) - A(y)| < |A(x) - Afo- n )| + |Afo- n ) - A(y)\ < 2 x S n (^ n ,e n ) 

• Otherwise and witout loss of generality, there exists j, k with k > j such that 
x e [£j, n ,€j+i,n] and y G [£k,n,£k+i,n\- Since \x - y\ < e n , it follows that |£ fc>n - 
£j+i,n| < £n- Then, 

|A(x) - A(y)| < |A(x) - A(&, n )| + |A(& >n ) - A& +1>n )| + |A& +1 , n ) - A(y)| 

In other words, for all x, y one may obtain 

\A(x) - A(y)\ < 3 x max e„). 

0<J<Pn 

Hence, S* < 3 x maxo<j< Pn S n (^j jn , e n ). Now, following the proof of Lemma 11, one 
may prove that there exists some positive constant 6 1 such that for n large enough and 
for all j = 0, . . . ,p n , 

P(5„(^,„,e„)> Vn) = 0(n"T). 
And in particular for 7 = 2, it comes 

P I max S n {£ j)n , e n ) > 2 r n I < (p n + 1) . max P (S n (£ jjn , e n ) > 6> 7 r n ) 

\0<J<Pn / J=0,...,p n 

= °(5)=°("" 3/2 )- 
whatever the value of ar POiP1 . This leads to the result by using Borel-Cantelli's Lemma. 
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6.5 Proof of Theorem 3 

Proof. We follow the proof of Theorem 2. Let p £ \po,Pi] and let e n = e n (a, r POiP1 ), 
then ^ 

P ~ f ^ ]) ~ (ftp) " f(P)) = A ^ + + C ^ 

where A(p),B(p) and C(p) are respectively defined by (48), (49) and (50). Similarly 
to the proof of Theorem 2, one may prove that sup po<p<pi A(p) = O a . s . { n ~ 1 )- Under 

Assumption A 5 (p , Pi), C(p) < (su PN < £n(Q)Tp) F"{x + g(p))) ■ Therefore, 

for n large enough, C(p) < A (sup po<p<pi (p) — £(p)^ • And from Lemma 12, this 



leads to 



sup C(p) = O a . s . {e n (a,T P0 , pi ) 2 ) 
po<p<pi 



In addition, using Lemma 13, one also has sup po <p< pi B(p) = O a . s . (r n (a, r POtPl )), which 
ends the proof. ■ 

6.6 Auxiliary Lemma for the proof of Theorem 4 
Lemma 14 Consider for < p < 1 the function h p (-), given by 

h P (t) = l{|t|<£ m (p)}(0 ~P, 



that is the function with g(-) = \ ■ \ . Then by denoting c* p the j-th Hermite 

coefficient of h p (-), we have for all j > 1 

4 P =% P +i = and c% = -2H 2j -i{q)<l>{q), (89) 
where q = i\ Y \(j)) = $- 1 (^). 

Proof. Since P (|Y| < q) = p and h p {-) is even, we have Cq p = c^ p - +l = 0, for all j > 1. 
Now, (27) implies 

4i = I h p (t)H 2j (t)(t>(t)dt = 2 x I" H 2j {t)4>(t)dt 
Jr Jo 

'^~^(t)] q = 2x[-H 2 ^ 1 (t)m q o 



= 2 x 



= -2H 2j -i(q)<i>(q). 
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Remark 10 Let g(-) = g(\ • \), where g(-) is a strictly increasing function on IR + ; then 
for all < p < 1, we have 

£|y|(p) = T X (t g (Y)(pj) • 

Consequently, the functions ^ 9(y) ( P )(") f or <7(") = I ■ \> <?(•) = I ■ \ a an d <?(•) = log | • | 
are strictly identical. And so, their Hermite decomposition is given by (89) and their 
Hermite rank is equal to 2. 

6.7 Proof of Theorem 4 
Proof, (i) Define 

M 

&„ = -^5 m log(l + <5f>)), (90) 

m=l 

where <5™ m (0) is given by (4). From (14), (15), and (25), we have almost surely 

M 

H a -H = 



M 

a m 

rn=l 



> log — i- + a x b n 



m=l 
M 



= E t \ (?> ^ \ Yam D ~ t\r\° CP. c )) (! + (!)) + « ( 91 ) 

^ <|y|- (P, c) V J 



and 



M 

plog 



-H=Y J B„ 

m=l 

M 

= Y, Bm (£(P> c ; lo § I 1 "*™ I) " ^log |^| (P, c)) + 6n (92) 

m=l 

Under Assumption Ae(f7), we have 

6 n = O (n"") . (93) 
Moreover, let i,j > 1, under Assumption A\(2v), we have, from Lemma 1 

V{Y am (i)Y am {i + j)) = p am (j) = O {\j\ 2H ~ 2 n • (94) 
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Then, for all m = 1, . . . , M and for all k = 1, . . . , K, from Lemma 10 and Remark 10, 
we obtain, that almost surely 

Z(Pk;\Y am \ a ) -£ m «(p fc ) = O a . a . (y n (2is-2H,T Pk )), 
£(p k ; log \Y am |) - 6og|Y|(Pfc) = O a . s . (y n {2v - 2H,r Pk )) , 

where the sequence y n (-, •) is defined by (52) with L(-) = 1. The result (37) is obtained 
by combining (91), (92) and (93). 

(ii) Let us apply Theorem 2 to the sequence g(Y am ), for some m = 1, . . . , M, with 
</(■) = | • |, g(-) = | • \ a and </(■) = log | • |. For all k = 1, . . . , K, we have almost surely 

£(p k ;\Y \)-Z\Y\{pk) = 7 77 7- u +0 . s .(r„) 

/|y|«(qy|(Pfc)J 

€ (p*; |v I ) - e,y|-(P*) = /|y|a(e|y|Q(pfc) ) + °~- (rn) 

f (ft; log |V" I) - 6 og ,y, (ft) = ^-^l^)^l ytt "l) + 0a . s . (r n ) , 

where, for the sake of simplicity, r n = r n (2z/ — 2H,T Pk ) defined by (35) and (36). Note 
that from Remark 10 T Pk = 2 for all k = 1, . . . , K. 

With some little computation, we can obtain, almost surely 

Hp k ;\Y am r)-Z\Y\«(Pk)=^\Y\{PkY~ l (UPk;\Y am \)-CiY\(Pk))+O a . s .(r n ), (95) 
and 

UPk^og\Y am \)-C loglY \(p k ) =Z\Y\(p k r 1 (UPk;\Y am \)-Z lYl (p k ))+O a . s .(r n ). (96) 

From (91), (92), (95), (96) and properties of Gaussian variables, the following results 
hold almost surely 

M K 

% C k 

II ■ J /, '! 

and 



R 



M K 



£iog _ H = £ (F( %; |y|) - Pfc ) + . s . (r n ) + O (6 n ) , (98) 

m=lfc=l qk< P^ k > 

where % and 7r£ are defined by (42). Denote by fc the following constant 

m > k 2q k ^q k ) k - 
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Since ir% = 1, (97) and (98) can be rewritten as 

H a -H = Z« + O a . s . (r n ) + O (b n ) (99) 
£iog _ H = z o + (rn) + Q {K) j (1Q0) 

where for a > 0, 

M K 

z n = EE^ 1*1) -Pk) ■ (ioi) 

m=l fc=l 

Thus, under Assumption A6(r?), we have, as n — > +oo, 

MSE(H a -H) = 0(B {(Z%) 2 )) + O (r n (2u - 2H, 2) 2 ) + O (n~ 2r >) , (102) 
MSE(H l °z -H) = 0(E ((Z°) 2 )) + O (r n (2v - 2H, 2) 2 ) + O (n" 2 ^) . (103) 

Now, 

M K 



E 



(ra 2 ) 4 E EE ( v e^ arai (^)) v 2 (^ am2 fe))) • 



n 

mi,m2=l fci,fe2=l ii,?2=l 



For &i, fc 2 = 1, . . . , if , mi, m 2 = 1, . . . , M and zi, i 2 = 1, • • • , i, we have from Lemma 14, 



E 



(V( ir ° mi («i))^( yttma (<=»))) = E E 



., (2ji)!(2j 2 )! 



X 



E(ii 2jl (y ami (n))ii 2j2 (y° m2 (z 2 ))) 



Co „' C. 



= E^yhp° mi ' om2 (^-n) 2j - (104) 

Under Assumption A\(2u), we have from Lemma 1, p ™ 1 ' " 12 (i) = O (\i\ 2H ~ 2u ) • Now, 
we leave the reader to check that, asn^+oo 

J_ £ ^ 1 ,a- ( . 2 _- i)2=0 | i ^ 02(217-2,) j =OK(2l/ _ 2 ^ )))j 
»i,»2=l \ |i|<n / 

where the sequence f n (-) is given by (39). Thus, we have, as n -► +oo, E((Z£) 2 ) = 
O (u n (2i/ - 2i2")), which leads to the result from (102) and (103). 

(in) Assume v > H + 1/4 and 77 > 1/2, then from (99) and (100), the following 
equivalences in distribution hold 

yfc(HZ-H)~yfc Z% and (h%* - H) ~ ( 105 ) 
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Now, decompose Z" = + T%, where 

M K Mt+l 

T « = 7^EE^ E v(^° m w) 

* m=l fc=l i=M-l 

and 

M K ( n ~\ 

^ 2 = ^EEc, fe U E vcr° ro «) , 

m=l fe=l I i=Jl«+l J 

Clearly, T^; converges to in probability, as n — ► +oo. Therefore, we have, as n — > +oo 

E Ga ( yal (<).-. ^ a "«)| (106) 

I i=M£+l J 
where G° is the function from R M to R defined for a > and ti, . . . , tu € R by: 

M K 

G a {t 1 ,...,t M ) = Y J Y. 9 ™*h qk {tm)- (107) 

m=l k=l 

Denote by Y a (i), the vector defined for i = Ml + 1, . . . , n by 

= (I- 1 (*),..., y- M (i)). 

We obviously have E(G a ("K a (i)) 2 ) < +oo. Since, for all k = 1,...,K, the functions 
h qk have Hermite rank r Pk , the function G a has Hermite rank 2 (see e.g. Arcones 
(1994) for the definition of the Hermite rank of multivariate functions). Moreover under 
Assumption Ai(2i/), we have from Lemma 1, as j — > +oo 

E (Y" mi (i)Y" m *(i + = p ami ^(j) = O (\j\ 2H ~ 2 ») € l\TL\ 

as soon as v > H + 1/4. Thus, from Theorem 4 of Arcones (1994), there exists cr 2 
(defined for a > 0) such that, as n — > +oo, the following convergence in distribution 
holds 

ZJ— «V(0,oS) 

with 

°a = E ^ ( G ° (^°(0) G° (W +%))). 

With previous notations, we have 

M K 

^=E E E c 1 , fcl c 2)fc2 E(v(^ 1 (0)^ 2 ^ am2 (^+^))) 

ieZ mi,m 2 =l fei,fe 2 =l 

M K hp H h Pk 2 

= E E E E^^^r^w 23 . (los) 

iSZ mi,m 2 =l fci,fc 2 =l j>r 
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From (89), we can see that formula (108) is equivalent to (41), which ends the proof 
from (105). ■ 

6.8 Proof of Corollary 5 

Proof. Equation (106) is still available for a sequence a n such that a n — > as n — ► +oo, 
that is 

{ n i=Me+i ' J 
From (107) and since 7r£ n — »■ 1, as n — >■ +oo, we have G Qn (-) — »■ G°(-). Therefore, the 
following equivalence in distribution holds, asn-> +oo 

which ends the proof. ■ 

6.9 Auxiliary Lemma for the proof of Theorem 6 

Lemma 15 Let < (3\ < (5?, < 1 and let Z = (Z\, . . . , Z n ) n random variables identi- 
cally distributed, such that sup ( g 1<p<( g 2 £z (p; Z) = O a , s , (1), then 

Proof. It is sufficient to notice that for i = 1, . . . , n 



zif3) - 1 r r t * & (p; z ) * = ( n_1 ) 

1 — PI — Pi J 0i 



j_ i + l 

/n — / n — 

i £ (P5 z ) dp < z (i),n < n I £ (p; Z) dp. 



which leads to 



n-[n0 2 ] n-[nf3 2 ] l 

— — - ; e (p; * < ^ (/3) < — nrr m L T " ^ z ) 

n - [n/? 2 J - [n/3ij 7 Mil n - [ra/?2j - [np\\ J Mil , i 



n n 



The end is omitted. ■ 

6.10 Proof of Theorem 6 
Proof, (i) From (20), (21), and (26), we have 
M B 



m=l 
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M _ 

£ £m log (|Y« ro |« (/3) /Pr (/3) ) + « x 6„ 

m=l 

M B. 



m=l \ a 



and 



M 

Jog, tm 



filog,tm_ H= ^ Bm e\ 
m=l 

M 

= £ s m (log |y« ro | (/3) - io?m (/3) ) + &n (no) 

m=l 

Let us notice that from Lemma 12, one can apply Lemma 15 for the vectors |i Aa ' Tl | c * 
and log \ Y am |. Then it comes 

if^ (/3) _|^) = 1 _l_ Pi f~" 2 (Uv, \r am r) - e, y |a ^ + o M . (n- 1 ) , 

log|^| (/3) -Io^m {/3) = 1 (p;log|y° m |) -6o g |y|(p)) dp + C.. K 1 ) 

Hence, from (94) and Lemma 12 and Remark 10 and under Assumption (Ae(rj)), we 
obtain 

fia,tm _ H = {yn{2H _ 2> 2)) + c + a s 

H los ' tm - H = O a . s . (y n (2H - 2, 2)) + O (n"") + ... (n" 1 ) , 

where the sequence y n (-, •) is defined by (52) with L(-) = 1. This leads to the result by 
noticing that n~ l = O (y n (2H - 2, 2)). 

(ii) By following the proof of Theorem 4 (ii) and from Theorem 3, we may obtain 
the following representation 



+ O . S . (r n ) + . s . (n' 1 ) + O , 

^F( q ;\Y" m \)-p 



H^tm_ H= Y j B m x 1 

^1 1 - P2 - A J/3! 



2#(< ? ) 

+ O a . s . (r n ) + . 8 . (n- 1 ) + O . 



dp 
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With such a representation, we observe that the result (ii) can be proved similarly to 
the one of Theorem 4. 

(Hi) By assuming that 77 > 1/2 and v > H + 1/4, one may obtain the asymptotic 
normality of # a ' tm and # lo g>* m by using the same tools as the one presented in the proof 
of Theorem 4 (in). Therefore, let us just explicit the asymptotic variance of estimators 
H a ' tm and H l ° s,tm . If v > H + 1/4 and rj > 1/2, then from previous representations and 
from 104 we obtain as in +00 

M n 1 1 



mi,m2=lii,t2=l (1^1 
■ 1-/32 .i-^ >i> a -l a -l 

hi hi j?[ ( 2 -?) ! tyKQimto) 

with % = for fc = 1, 2. Due to (89) and since \Y^ {(3) = J^ 2 q a dp, 

this variance converges towards <Jq j m given by (46), as n — > +00. 

We leave the reader to check that the asymptotic variance of y/n ^H log ' tm — H^j is 
given by al tm . ■ 
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Figure 1: Two examples for the sample paths of non-contamintcd (top) and contaminated 
processes with variance function v(-) = \ ■ \ 2H (left), respectively v(-) = 1 — cxp(— | • \ 2H ) (right), 
see (47). 
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Figure 2: Left: cr^ tm in terms of f3; Right: o\ for estimators based on a single quantile 
in terms of p. Three values of the parameter H are considered: 0.3 (top), 0.5 (middle), 
0.8 (bottom). The parameter M is fixed to M = 5. The constant line corresponds to 
the asymptotic variance of the Whittle's estimator. 
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Non-contaminated sample paths 



Estimators „(•) = | • \ 2H v(-) = l - ex P (-| ■ \ 2H ) 

p = 1/2, c = 1 (median) 0.796 (0.042) 0.801 (0.042) 

p = 0.9, c = 1 0.797 (0.035) 0.798 (0.036) 

p = (1/4, 3/4),c = (1/2, 1/2), g(-) = | • | 2 0.795 (0.036) 0.800 (0.037) 

10%-trimmed mean, g(-) = | • | 2 0.797 (0.03) 0.799 (0.034) 

Quadratic variations method 0.802 (0.032) 0.798 (0.032) 

Whittle estimator 0.805 (0.024) 0.806 (0.024) 



Contaminated sample paths 



Estimators «(•) = | • \ 2H v(-) = l - ex P (-| ■ \ 2H ) 

p = 1/2, c = 1 (median) 0.798 (0.047) 0.803 (0.045) 

p = 0.9, c = 1 0.793 (0.033) 0.789 (0.032) 

p = (1/4, 3/4),c = (1/2, 1/2), g(-) = | • | 2 0.797 (0.040) 0.796 (0.037) 

10%-trimmed mean, g(-) = | • | 2 0.792 (0.037) 0.797 (0.033) 

Quadratic variations method 0.329 (0.162) 0.353 (0.149) 

Whittle estimator 0.519 (0.106) 0.510 (0.100) 



Table 1: Mean and standard deviations for n = 1000 and H — 0.8 using 500 Monte Carlo 
simulations of sample paths of processes with variance function v(-) = \ ■ \ 2H , respectively v(-) = 
1 — exp(— | • \ 2H ) (first table ) and contaminated versions (second table), see (47). 



